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In this paper we introduce a perturbatively super-renormalizable and unitary theory of quantum 
gravity in any dimension D. In four dimensions the theory is an extension of the Stelle higher 
derivative gravity that involves an infinite number of derivative terms characterized by two entire 
functions, a.k.a. "form factors". In dimension D we preserve two entire functions and we implement 
a finite number of local operators required by the quantum consistency of the theory. The main 
reason to introduce the entire functions is to avoid ghosts (states of negative norm) like the one 
in the four-dimensional Stelle's theory. The new theory is indeed ghost-free since the two entire 
functions have the property to generalize the Einstein-Hilbert action without introducing new poles 
in the propagator. By expanding the form factors to the lowest order in a mass scale we introduce, 
the local high derivative theory is recovered. Any truncation of the entire functions gives rise to the 
unitarity violation and it is only by keeping all the infinite series that we overcome similar issues. 
The theory is renormalizable at one loop and finite from two loops upward. More precisely, the 
theory turns out to be super-renormalizable because the covariant counter-terms have less derivatives 
then the classical action and the coefficients of the terms with more derivatives do not need any 
kind of infinity renormalization. In this paper we essentially study three classes of form factors, 
systematically showing the tree-level unitarity. 

We prove that the gravitation potential is regular in r = for all the choices of form factors 
compatible with renormalizability and unitarity. We also include Black hole spherical symmetric 
solutions omitting higher curvature corrections to the equation of motions. For two out of three 
form factors the solutions are regular and the classical singularity is replaced by a "de Sitter-like 
core" in r = 0. 

For one particular choice of the form factors, we prove that the _D-dimensional "Newtonian cosmol- 
ogy" is singularity-free and the Universe spontaneously follows a de Sitter evolution at the "Planck 
scale" for any matter content (either dust or radiation). 

We conclude the article providing an extensive analysis of the spectral dimension for any D and 
for the three classes of theories. In the ultraviolet regime the spectral dimension takes on different 
values for the three cases: less than or equal to "1" for the first case, "0" for the second one and 
"2" for the third one. Once the class of theories compatible with renormalizability and unitarity is 
defined, the spectral dimension has the same short distance "critical value" or "accumulation point" 
for any value of the topological dimension D. 



In 1916 Albert Einstein revolutionized the way physi- 
cists thought about gravity with the theory of "general 
relativity" . This theory works quite well at the classical 
level, but at the theoretical level one of the biggest prob- 
lems is to find a theory that is able to reconcile general 
relativity and quantum mechanics. There are many rea- 
sons to believe that gravity has to be quantum, some of 
which are: the quantum nature of matter in the right- 
hand side of the Einstein equations, the singularities ap- 
pearing in classical solutions of general relativity, and so 
on. 

Let us recall here the main hypothesis underlying gen- 
eral relarivity. The grounding principles of "general rel- 
ativity" are: (i) gravity is no longer a force as in the 
Newton's theory, but it is the "curvature of the space- 
time" , (ii) the symmetry principle underlying the gravi- 
tational theory is the "general covariance" or "invariance 
under general coordinate transformations" , (iii) the "en- 
ergy momentum tensor" for any kind of matter has to be 
covariantly conserved, (iv) the dynamics should be de- 
scribed by "second order" differential equations. Assum- 
ing these fundamental requirements, we get a "unique" 
theory for the dynamics of gravity, namely the Einstein 
equations. The theory can be also formulated starting 



from an action principle by Einstein and Hilbert. The 
action principle we are going to introduce in this paper 
makes gravity compatible with quantum mechanics, as 
it is the result of a synthesis of minimal requirements 
beyond the classical Einstein-Hilbert action. 

Let us axiomatically itemize these requirements one by 
one: 

(i) . gravity is still retained as curvature of the space- 

time and the underlying symmetry principle re- 
mains "general covariance" ; 

(ii) . Einstein-Hilbert action is no longer the correct one, 

but it should still be a good approximation of the 
theory at an energy scale much smaller than the 
Planck mass or any other invariant scale; 

(iii) . solutions of the classical equations of motion must 

be singularity-free; 

(iv) . the theory has to be perturbatively renormalizable 

or finite at quantum level. In other words, we as- 
sume "axiomatically" the validity of the "perturba- 
tive theory" . We claim that any mathematical the- 
ory which prides itself to describe the universe must 
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be perturbatively consistent. This is empiricahy 
true for ah the other fundamental forces: Weak, 
Strong and Electromagnetic. If for a general sys- 
tem we find that this is not the case, then we have 
to change the theory or "dualize" it [iHll instead 
of trying to solve it at a non-perturbative level; 

(v) . the spacetime spectral-dimension should decrease 

with the energy to obtain a complete quantum 
gravity theory in the ultraviolet regime (this hy- 
pothesis is strongly related to the previous one 
(iv)). This property seems to be of "universal na- 
ture" because it is common to many approaches to 
quantum gravity The Stelle's theory [l^, 

the Crane-Smolin theory "asimptotically safe 
quantum gravity" |19l - [25| , "causal dynamical trian- 
gulations" [l^ , "loop quantum gravity" Q and 
"string theory" already manifest this property with 
an high energy spectral dimension rf^ = 2 or less; 

(vi) . the theory has to be unitary, with no other pole 

except the graviton one in the propagator; 

(vii) . spacetime is a single continuum of space and time 

and in particular the Lorentz invariance is not bro- 
ken. This requirement is supported by recent ob- 
servations. 

The main hypothesis we abandon with respect to the 
classical theory is the absence of higher derivative terms 
in the action. As we are going to show we admit indeed 
an infinite tower of operators defined through an "entire 
function" of the D'Alembertian differential operator. 

This work is inspired by papers in four dimensions 
about a "nonlocal extension" of gauge theories and grav- 
ity published in the Nineties [23, [33 - l39j . For exam- 
ple, in [33l - l39| the authors considered a modification of 
the Feynman rules, where the coupling constants {gi for 
electro-weak interactions and Gn for gravity) become 
functions of the momentum p. They proved the gauge 
invariance at all orders in gauge theory but only up to 
the second order in gravity. For particular choices of 
gi(j)) or Gn{p), the propagators do not show any other 
pole besides those of the standard particle content of the 
theory, which makes the theory unitary. On the other 
hand the theory is also finite if the coupling constants 
go sufficiently fast to zero in the ultra-violet limit. On 
the basis of these conclusions, the problem with gravity 
remains to find a covariant action that self-contains the 
properties mentioned above: finiteness and/or renormal- 
izability and unitarity. 

We now introduce our theory starting from the per- 
turbativc ZJ-dimensional "non-renormalizable" Einstein 
gravity, and recalling the four dimensional Stelle's 
quantum gravity, which serves as an example of power- 
counting renormalizable (but not unitary) theory. 

First of all, let us explain why it is important to study 
quantum gravity in Z3-dimensions. There are at least 
five reasons to look for a super-renormalizable theory of 
gravity with extra dimensions, (i) The first reason is 



that D > 4 eliminates "soft-graviton or infra-red diver- 
gences" at quantum level [i^l; (ii) the second one is that 
only in £) ^ 7 there exists a well defined "total cross sec- 
tion" [i^l; (iii) the third one is to obtain a well defined 
"Kaluza-Klein grand-unification" ; (iv) the forth one is 
the possibility to have a well defined completion of the 11- 
dimcnsional supcrgravity as a candidate for "M-Theory" ; 
(v) the last reason to study gravity in any spacetime di- 
mension is related to the universality of the quantization 
procedure independently from the number of dimensions 
"Z3" . In synthesis, we do not want to tune the spacetime 
dimension to a particular value to make the theory con- 
sistent at quantum level; instead, the theory "should be" 
well defined for any value of D. 

In this paper we use the signature (H — ... — ). The 



curvature tensor is defined by i?^ 



'MS 



-dsn. 



the 



Ricci tensor by i?^^ ~ R'^ua i ^^'^ the curvature scalar by 
R = g^^'^ Rfj^n and g^j/, where the metric tensor. 

Perturbative quantum gravity is the quantum theory 
of a spin-two particle on a fixed (conventionally flat) 
background. Starting from the _D-dimensional Einstein- 
Hilbcrt Lagrangian 



C 



EH 



g\2K-'R, 



(1) 



(k^ = 327rG7v) we split the metric in a background part 
and in a fluctuation 



Kh,, 



(2) 



Using ([2]), we then expand the action in powers of the 
graviton fluctuation h^^^, around the fixed background g"^ 
that we assume to be the flat Minkowski metric ?7^i/. Re- 
grettably, the quantum theory is divergent at two loops 
in D = 4, so we should introduce a counter-term propor- 
tional to the Riemann tensor at the third power 



£^3 cx G 



N ■ 



(3) 



In the general D dimensional Einstcin-Hilbert theory the 
superficial degree of divergence of a Feynman diagram is 
S ^ L D + 2V — 21, where L is the number of loops, V is 
the number of vertices and / the number of internal lines 
in the graph. Substituting in S the topological relation 
L = 1 + I — V. wc obtain 



S = 2 + {D -2)L for C 



EH ■ 



(4) 



In _D > 2 the superficial degree of divergence increases 
with the number of loops. Therefore, we are forced to 
introduce an infinite tower of higher derivative counter- 
terms and an infinite number of coupling constants, thus 
making the theory not predictive. 

Schematically, we can relate the number L of loop di- 
vergences to the counter terms wc have to introduce to 
regularize the theory. In short 



5 = / d^ 



XV \9 



4-CK3 



(5) 



n+2m=2+{D-2}L 
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where "n" and "m" are integer numbers, amn are cou- 
pling constants and 1/e is the cutoff in dimensional reg- 
ularization. 

Now we want to comment further about the above 
statement on the meaning of renormalizability. Scholars 
usually claim that a theory loses its predictability when 
its action consists of an infinity number of operators. 
Consequently, they believe that the theory can be defined 
only through an infinite number of measures. Such state- 
ment is highly questionable because we should always 
add all the possible operators to a Lagrangian and/or 
an Hamiltonian describing a physical system. The most 
important step we need to take is to assess whether the 
physical measurable quantities are affected or not by the 
above-mentioned operators. In other words it is essen- 
tial to establish if such operators are "relevant o irrele- 
vant" to the physical observable. If we can assume that 
no physical quantity is susceptible to such operators, we 
can then empirically infer that the coupling constants 
equal "zero", as long as the theory satisfies the unitar- 
ity requirement. Let us provide an example. Suppose to 
add other hcrmitian operators to the standard Hamilto- 
nian of the hydrogen atom. We can "invent" an infinite 
number of such operators but only few of them will be 
relevant, like for example the relativistic corrections or 
the Lamb Shift, while all the other ones are irrelevant. 
In this case what is really significant is the precision level 
in the measure of the energy spectrum, compatibly with 
a zero value for the other couplings, as long as unitarity 
is satisfied. 

A first revolution in quantum gravity in 4Z) was intro- 
duced by Stelle [itI . |4]| with his higher derivative theory 



2k-^R 



(6) 



If we calculate the upper bound to the superficial degree 
of freedom for this theory in _D-dimensions we find 

5 = DL- 41 + W = D - {D - 4){V - I) , (7) 

so that (5 = 4 in £> = 4. The theory is then renormal- 
izable, since all the divergences can be absorbed in the 
operators already present in the Lagrangian 1^. Un- 
fortunately however, the propagator contains a physical 
ghost (state of negative norm) that represents a violation 
of unitarity. Probability, as described by the scattering 
S'-matrix, is no longer preserved. Similarly, the classical 
theory is destabilized, since the dynamics can drive the 
system to become arbitrarily excited, and the Hamilto- 
nian constraint is unbounded from below. On this basis, 
we can generalize the Stelle theory to a i'-dimensionl 
renormalizable one. In short, the Lagrangian with at 
most X derivatives of the metric is 



^D-Rcn — O.1R + a2R^ + b2R^^ + 

... + a^R""^ + b^R'^l^ + cxR^l. + dxRU'i-^R . 



(8) 



In the second line, the dots on the left imply a finite num- 
ber of extra terms with fewer derivatives of the metric 



tensor, and the dots on the right indicate a finite num- 
ber of operators with the same number of derivatives but 
higher powers of the curvature {0{R'^^'^/^~'^R)). 

In this theory, the power counting tells us that the 
maximal superficial degree of divergence of a Feynmann 
graph is 



6 = D - [D - X){V - I). 



(9) 



For X = D the theory is renormalizable since the max- 
imal divergence \s 5 ^ D and all the infinities can be 
absorbed in the operators already present in the action 

(ED. 

The general action of "derivative order A^" can be 
found combining curvature tensors and covariant deriva- 
tives of the curvature tensor. In short the action reads 
as follows [42} . 

S^Y. "2„A''~'" / d''x^\02n{dpg^,) + 5nl , (10) 

n=0 '' 

where A is a mass scale in our fundamental theory, 
02n{dpg^u) denotes the general covariant scalar term 
containing "2n" derivatives of the metric .g^,y , while S'nl 
is a nonlocal action term that we are going to set later 
[isl - fisj . The maximal number of derivatives in the local 
part of the action is 2A -I- 4. We can classify the local 
terms in the following way. 



A^ = 1 : 5i = 5*0 + c^i^R^„ 



JO] 02 



„(0) 



C2 'WR..yR.. 



N = 2 : S2 = So + Si+ c-pRt 

+c^2^R..yR..yR.. 



.(3)Y72i 



C3^V^i?...V^i?.. 



N-1,N>0 



N ^ N 



3JV 



i=0 



J AT) jyN-lr 



-C2 'i?;:-'vi?...vi?. 



-cWi?...n^i?... . (11) 



In the local theory ([8]), renormalizability requires X = D, 
so that the relation between the spacetime dimension and 
the derivative order is 2 A + 4 = D. To avoid fractional 
powers of the D'Alembertian operator, we take 2A^-|-4 = 
Dodd -l- 1 in odd dimensions and 2A^ + 4 = Devon in even 
dimensions. 

In this paper, we are focused on the renormalizability 
and unitarity of the theory, so the main quantity to cal- 
culate is the graviton propagator. Although the action 
is complicated, we only need to consider the quadratic 
operators in the curvature to get the "two points func- 
tion". Given the general structure ([TU]). for A^ > and 
n > 2 contributions to the propagator come only from 
the following operators. 
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However, using the Bianchi and Ricci identities one can 
reduce the terms hsted above from three to two (with 
total 2n derivatives) pl |. 



(13) 



= m^,^^^-'^R^"' - Ra"-^R + 0{R^) + V^l^'^ , 

where V^fi^ and V^51'^ arc total divergence terms. Ap- 
plying ([T3| to (fT2|). for n > 2 we discard the third term 
and we keep the first two. 

We now have to define the "non-local" action term 
in ()10p . As we are going to show, both supcr- 
rcnormalizability and unitarity require the following two 
non-local operators. 



R^,h2{-aA)R^\ Rho{-nA)R. 



(14) 



The full action, focusing mainly on the non-local terms 
and the quadratic part in the curvature, reads 



S = J d^'x^l 

N 



2k-^R + X 



-Y,[anR (-Oa)" R + bn R^, (-Da)" i?^" 

n=Q 

-Rhni-OA) R + i?^. h2{-aA) R^''' 



+ . 



.0(ii'3). 



R 



N+2 



(15) 



Finite number of terms 



The last line is a collection of local terms that are renor- 
malized at quantum level. In the action, the couplings 
and the non-local functions have the following dimen- 
sions: [a„] = [bn] = M"-^, [k^] = AP-^, [A] = A/^, 
[h2] = [ho] = 

At this point, we are ready to expand the Lagrangian 
at the second order in the graviton fiuctuation. Splitting 
the spacctimc metric in the fiat Minkowski background 
and the fiuctuation /i^i^ as defined in ([2]), we get [4^ 

Clin = ^^[h^^'^Dh^, + Al + {A, - 0^,)2] 



+2k^ 00) «(□)(< -□(/.) 



(16) 



where = /i^^^,, (j) ^ h (the trace of h^^y), F^i, = A^^^y — 
Ajj^f^ and the functionals of the D'Alembertian operator 
/3(n),a(n) are defined by 

N 

ri=0 
N 

/3(n)/2:=^&„(-nA)" + /i2(-nA). (i7) 



The d'Alembertian operator in ^TE\\ and pT)) must be 
conceived on the fiat spacetime. The linearized La- 
grangian (fT6|) is invariant under infinitesimal coordinate 
transformations ^ + k^^{x), where ^''(2;) is an in- 
finitesimal vector field of dimensions [^(a;)] = M'^^^^^^^. 
Under this transformation, the graviton field turns into 



(18) 



The presence of the local gauge simmetry ([T8| calls for 
the addition of a gauge-fixing term to the linearized La- 
grangian ()16p . Hence, we choose the following fairly gen- 
eral gauge-fixing operator 

Cgf = Ai(A, - A0,,)c^i(-nA)(A'^ - Xcb^n 
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-(Af; - An0)/?(n)c.2(-nA)(A;; - ad^) 



/3(n)c^2(-nA) 



(19) 



where a;i(— Da) are three weight functionals [13, S^- 
the harmonic gauge A = A2 = A3 = and Ai = 1/^. The 
linearized gauge-fixed Lagrangian reads 



^lin + ^GF — -jh^^ O ^u^pah'"^ , 



(20) 



where the operator O is made of two terms, one coming 
from the linearized Lagrangian (|16p and the other from 
the gauge-fixing tcrm (|T9)) . Inverting the operator O (43| 
we find the two point function in the harmonic gauge 

(5^V = 0)> 



p(0) 

2k' 



D/3{k^)K^/i+{D-l)a{k^)K2 
2 



(21) 



The tensorial indexes for the operator O ^ and the pro- 
jectors p(o),p(2),p(i),p(o) have been omitted and the 
functions a{k'^) and I3{k^) arc achieved by replacing 
— □ — > k^ in the definitions (|17p . The projectors are 
defined by [H [gI 

P^wpai^) " 2 (^MP'^I^'T + Qfia^^vp + dupLOp^cr + Qua^^fip) , 
Pjiupai^) = -jy—^Opudpa , P^xupai^) ^ ^ pv^ pa ^ 



^ pupa " pv^pa + '-^pu"pa i 

p. kf^ki, 

^ pf Vpi' 7 ^ pu 



knkii 



(22) 



Looking at the last two gauge invariant terms in (|2ip . 
we deem convenient to introduce the following defini- 
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tions, 



n=0 



:h2iz).. 



(23) 



D-2 



ho{z) 



D-2 K^A^D 



n=0 
N 



-K^A^{D - l)z\^^anz'' + ho{2 



where again z = —Da. Through the above definitions 
(j23[) , the gauge invariant part of the propagator greatly 
simplifies to 



1 /P(2) 



p(0) 



h2 [D - 2)ho 



(24) 



We clarify now the incompatibility of the unitarity 
with a polynomial choice of either /io,/i2 or, equally, of 
the two functions a(n), 13(D). If we assume for a moment 
a{D) and /?(□) to be polynomial, then each of the two 
fractions in the propagator takes the following form, 



1 



fc2(l+p„(fc2)) k 



where p„(x) is a polynomial of degree n. In (|25p we 
used the factorization theorem for polynomials and the 
partial fraction decomposition psl ]. When multiplying 
the left and right side of by and considering the 
ultraviolet limit fc^ — >■ +oo, we find that at least one 
of the coefficients Ci is negative. Therefore the theory 
contains at least a ghost in the spectrum. The conclusion 
is that h2,ho cannot be polynomial. 

Once established that /12 and hg are not polynomial 
functions, we demand the following general properties 
for the transcendental entire functions hi{z) [i = 0,2) 
and/or h,{z) [i = 0,2) [H: 

(i) hi{z) [i = 0, 2) is real and positive on the real axis 
and it has no zeroes on the whole complex plane 
l^l < +00. This requirement implies that there are 
no gauge-invariant poles other than the transverse 
massless physical graviton pole. 

(ii) 1^.1(2)1 has the same asymptotic behavior along the 
real axis at ±00. 

(iii) There exists > such that 

lim \h,{z)\^\zy+^, 

i-l-00 



-1^D/2 for i? = i?cvon, 

7 ^ (i? - l)/2 for D = Dodd 



(26) 



for the argument of z in the following conical re- 
gions 

C = {z\ - e < argz < -1-9 , tt - 6 < argz < tt + 9}, 
for < e < tt/2. 



This condition is necessary in order to achieve the 
(supe-)renormalizability of the theory that we are 
going to show here below. The necessary asymp- 
totic behavior is imposed not only on the real axis, 
(ii) but also in the conic regions that surround it. 
In an Euclidean spacetime, the condition (ii) is not 
strictly necessary if (iii) applies. 

Let us then examine the ultraviolet behavior of the quan- 
tum theory. According the property (iii) in the high en- 
ergy regime, the propagator in the momentum space goes 
as 

0-i(fc)^l/fc27+2^+4 for large ^ 

(see [231 mi)). However, the n-graviton interaction 
has the same leading scaling of the kinetic term, since it 
can be written in the following schematic way, 

£(") ^ /i"n^/i h,{-nA) a^h 

/i" D,jh (□,, + h"' dhdy+^ D^h, (27) 

where the indexes for the graviton fluctuation h^^, are 
omitted and /ii(— Da) is the entire function defined by 
the properties (i)-(iii). From ([?7)) . the superficial degree 
of divergence in a spacetime of "even" dimension is 

<5cvc„ = i^cvcni - (27 + 2N + 4)1 + (27 + 2N + 4)V 

even 

)/+(27 + L> 

even 

= Z?cvc„-27(L-l). (28) 

On the other hand, in a spacetime of "odd" dimension 
the upper limit to the degree of divergence is 



Sodd = Dodd-{2-/ + l){L-l). 



(29) 



In (|28p and (|29| we used again the topological relation be- 
tween vcrtcxes V, internal lines / and number of loops L: 
I = V + L-1. Thus, if7>i?evcn/2or7> (I?odd-l)/2, 
only 1-loop divergences exist and the theory is super- 
renormalizable^. Only a finite number of constants are 
renormalized in the action (|15p . i.e. k. A, a„, 6„ and the 
finite number of couplings that multiply the operators in 
the last line. The renormalized action reads 



2Z^,K~^ R + Z-.X 



5= / d^x^\g\ 

N 



- J2 [Za,. a„ R (-Da)" R + Zk„ 6„ i?^. (-Da)" i?^' 

-i? hoi-OA) R + R^, /i2(-nA) i?^" 
Z a,cf'i?3 + + c[''^R^+' 



(30) 



^ A "local" supcr-renormalizable quantum gravity with a large 
number of metric derivatives was for the first time introduced in 

El. 
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All the couplings in (|30|) must be understood as renor- 
malized at an energy scale ^. Contrarily, the functions 
hi(z) are not renormalized. To understand this point 
thoroughly, we can write the generic entire functions as 
a series, i.e. hi{z) = '^^^a.rz'^. Because of the su- 
perficial degrees of divergence ([28| and (|29l). there are 
no counterterms that rcnormalizc for r > N. As a 
matter of fact, the couplings in the second line of ([501 al- 
ready incorporate the renormalizations of the coefficients 
Ur for r ^ N. Therefore, the non-trivial dependence of 
the entire functions hi{z) on their argument is preserved 
at quantum level. 

Imposing the conditions (i)-(iii) we have the freedom 
to choose the following form for the functions hi, 



holz) 



hoiz) 



V^(z)-1-1 + k^A^zELo«- 

K^A^ z 



, (31) 



for general parameters a„ and 6„. Here V{z)~^ = e^^'^^ 
and H[z) is an entire function that exhibits logarithmic 
asymptotic behavior in the conical region C. The form 
factor exp H{z) has no zeros in the entire complex plane 
for |z| < +00. Furthermore, the non-locality in the action 
is actually a "soft" form of non locality, because a Taylor 
expansion of hi{z) eliminates the denominator Da at any 
energy scale. 

The entire function H{z), which is compatible with the 
property (iii), can be defined as 



H(z) 



(32) 



where P'y+Ar+i(z) and C(z) must satisfy the following re- 
quirements: 

a. p^+N+iiz) is a real polynomial of degree 7 -I- + 1 
with p-f+N+i{0) = 0, 

b. C(^) is an entire and real function on the real axis 
with C(0) = 1, 

c. \C{z)\ — > for \z\ — > 00 in the conical region C 
defined in (iii). 

Let us investigate the unitarity of the theory. We as- 
sume that the theory is renormalized at some scale fiQ. 
If we set 



an = an(Mo) , bn = fon(Mo), 



(33) 



the bare propagator does not possess other gauge- 
invariant pole than the physical graviton one. Since the 
energy scale fig is taken as the renormalization point, 
we get h2 = ho = V{z)~^ = expiJ(z). Thus, only the 
physical massless spin-2 graviton pole occurs in the bare 
propagator and reads 



fc2 



p(2) _ 



p(0) 



D 



(34) 



The momentum or energy scale at which the relation 
between the quantity computed and the quantity mea- 
sured is identified is called the subtraction point and is 
indicated usually by "/i" [s^l. The subtraction point is 
arbitrary and unphysical, so the final answers do not have 
to depend on the subtraction scale fi. Therefore, the 
derivative d/d^^ of physical quantities has to be zero. In 
our case, if we choose another renormalization scale fi, 
then the bare propagator acquires poles. However, these 
poles cancel in the dressed physical propagator because 
the shift in the bare part is cancelled by a corresponding 
shift in the self energy. The renormalized action ([50)1 will 
produce finite Green's functions to whatever order in the 
coupling constants we have renormalized the theory to. 
For example, the 2-point Green's function at the first or- 
der in the couplings a„, 6„ can be schematically written 
as 



G2R = Vik^A') {e - EH(fc')), 



(35) 



where the renormalization prescription requires that S 
satisfies (on shell) 



T.r{Q) = and 



= 0. 



(36) 



We did not consider the tensorial structure and the lon- 
gitudinal components because they project away when 
attached to a conserved energy tensor. 

The subtraction point is arbitrary and therefore we can 
take the renormalization prescription off shell to = jj? . 
The couplings we wish to renormalize must be depen- 
dent on the chosen subtraction point, anifJ-) and bn{fJ-), 
in such a way that the experimentally measured cou- 
plings do not vary on shell. The renormalized Green's 
function G^^ at fj? must produce the same Green func- 
tion when expressed in terms of bare quantities. Conse- 
quently, the scalings Za,^ and Zh^ must also depend on 
/i^ . The Green's function written in terms of bare quan- 
tities can not depend on /i^, but since p^ is arbitrary, 
the renormalized Green's function must not either. This 
fact. 



d^I^ 



0, 



(37) 



is well known as the renormalization group invariance. 

When h2{z) = hQ{z) ~ 0, the action in reduces 
to ^ and fails to be unitary. Unitarity can be achieved 
only if 



in dll). 



which corresponds to the D-dimensional Einstcin-Hilbcrt 
non rcnormalizable action. 

An explicit example of exp H{z) that satisfies the prop- 
erties (i)-(iii) can be easily constructed. There are of 
course many ways to choose ((z), but we focus here on 
the exponential choice ^(z) ~ cxp(— z^), which satisfies 
requirement c. outlined for p2l) in a conical region C 
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with = 7r/4. The entire function H{z) is the rcsuh of 
the integral defined in ([5^ 

H{z) ^^[lE+r (0,p2^jv+i(^))] + log (p-t+N+iiz)) , 
Re(p^^^+i(z)) > 0, (38) 

where je ~ 0.577216 is the Eulcr's constant and 
r(a, z) = J^°° t°'~^e~^dt is the incomplete gamma func- 
tion. If we choose p^+N+i{z) = H{z) simplifies 
to: 



H{Z) = \[lE+T (0, ^ log(,7+A^+l) , 



(39) 



Another equivalent expression for the entire function 
H{z) is given by the following series 



+ 00 



H{z)^Y.^-iy 



,2n 



n=l 



2nn\ 



Rc(p^+jv+iW)>0. (40) 

For p^+M+i{z) = zT+^+i the 9 angle, which defines the 
cone C, is = 7r/(47 + AN + 4). According to the above 
expression (j40p we find the following behavior near z = Q 
for the particular choice p~f+N+i{z) = z''^^^^ , 



H{z) 



(41) 



We now present a systematic study of the tree-level 
unitarity [ijl. A general theory is well defined if 
"tachyons" and "ghosts" are absent, in which case the 
corresponding propagator has only first poles at fc^ — 
AP = with real masses (no tachyons) and with positive 
residues (no ghosts) . Therefore, to test the tree- level uni- 
tarity of a multidimensional super-renormalizable higher 
derivative gravity we couple the propagator to external 
conserved stress-energy tensors, T'"^, and we examine the 
amplitude at the pole. When we introduce a general 
source, the linearized action including the gauge-fixing 
reads 



hT 



The transition amplitude in momentum space is 



(42) 



(43) 



where g is an effective coupling constant. Here, only the 
projectors P^^^ and will give a non zero contribution 
to the amplitude since the energy tensor is conserved. 
To make the analysis more explicit, we can expand the 



sources using the following set of independent vectors in 
the momentum space [43j . 



— (^k ^ k'j ^ — i^k ^ k^ ^ 



(44) 



where ei are unit vectors orthogonal to each other and to 
k. The symmetric stress-energy tensor reads 



T"" ^ akt'k'' + hk^k" 



(45) 



The conditions fc^T'"' = and k^k^T''" = place con- 
strains on the coefficients a, b, d, e\ /* fisj . 

Introducing the spin-projectors ((22|) and the conserva- 
tion of the stress-energy tensor k^T'^^ = in (|43p . the 
amplitude results 



D 



Clearly, there is only the graviton pole in fc^ 
residue at fc^ = is 



(46) 
and the 



Res (A) 1,2 „ 



D-2 



\k-=o- (47) 



For D > 3 the result above result tells us that 
Res {A) |;,2_Q > 0, which means that the theory is uni- 
tary. Instead, for D — 3 the graviton is not a dynamical 
degree of freedom and the amplitude is zero. 

A first example of this quantum transition is the in- 
teraction of two static point particles. In this case, 
T^^ = diag(p, 0, 0, 0) with p — M 6{x) and the amplitude 
simplifies to 



A 



2 2 

9 P 



D-3 
D-2 



(48) 



which is positive in > 3 and zero for D = 3 since, 
again, there are no local degrees of freedom in D = 3. 

A second example we want to consider is the light 
bending in the multidimensional nonlocal gravity. We 
consider a static source and a light ray. The amplitude 
for this process is 



(49) 



where T'^" is the above energy tensor for the static parti- 
cle and is the traceless electro-magnetic energy ten- 
sor. Using the projectors defined in ((22l) and the propa- 
gator we obtain 



A = g'r'"T[ 



2rp00rp00 

EM 



(50) 



We see that, at low energy (A — > +oo), the amplitude 
(|50p is precisely the amplitude for the interaction between 
a static source and a light ray in £)-dimcnsional linearized 
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To address the problem of classical singularities men- 
tioned at the beginning of the paper, we can start out by 
calculating the gravitational potential. Given the mod- 
ified propagator p4[) . the graviton solution of the equa- 
tions of motion resulting from the Lagrangian (|42p with 
g = k/2, is 



= - I d^-r' I ^ ^ ^tk(x-x') 



D 



■T 



(52) 



fc2 V"^"- 
For a static source with energy tensor 

r,^ = diag(Af<5^-i(f),0,...,0), (53) 
the solution of ([5^ satisfying spherical symmetry reads 



-E,. 



(27r) 



D-l 



kM E, 



2 (27r) 



D-2 „D-3 



fc2 



D-1 p2 
2 '^T 



(54) 



where oFi{a; z) = oFi{a; z)/T{a) is the regularized hy- 
pergeometric confluent function. In (|54p . we also have 
introduced the variable p = \k\r and the matrix 



E, 



diag 



1 



D-2, 



D-2' (L>-2) 



1 



D-2 



Einstein's gravity. On the other hand, at high energy the 
light bending is much smaller in the nonlocal theory than 
in the Einstein's one. 

The proposed theory is not unique, but all the freedom 
present in the action can be read in the entire function 
V{z) or H{z) mils,!!!-!!!. The expression "form fac- 
tor" for the function V{z) used throughout the paper is 
not accidental. Indeed, It may be read in analogy with 
the form factors present in the interaction between a pho- 
ton and a nucleon. Most importantly, in view of this, it 
can also be measured experimentally. Specifically, the 
quantities and ([SO]) , at least in principle, can be 
measured experimentally to fit V{z). The four graviton 
amplitude is another measurable quantity of the form 
factor V{z) and it reads 

^4-grav. = A-grav. (s, i, U] V{s, t, u); 61,2,3,4) , (51) $(r) = 

where £1,2,3,4 are the four gravitons polarizations and 
s,t,u the Mandelstam variables. Since V{z) has to be = 
an entire function, we can falsify the theory by compar- 
ing the experimental four- gravitons amplitude with (|5ip . 



For r — J- 0, the entire function H{z) w log 2:'^+^+i and 
the solution (|54p is approximated by 



D-3 



2 (271)^-2^^''^ 

D-2 p2 
4 



X |dpp^-(2A^+4K2.-2^^^(^£_2._^^ 



(55) 



The solution ([55]) is clearly regular near r sa since the 
exponent of the radial coordinate is always positive in 
any dimension D. 

The gravitational potential is related to the /iqo com- 
ponent of the graviton field by $ = k/ioo/2. Then, using 
(|54|) we get 



4 £1-2 7 (27r)^-i 
GmM T(^] 



rD-3 



D-2 



D-3 

n 2 



■)4-Z3 



dp ■ 



\ 2 

For example, in D — 4, (|56p simplifies to 
GnM 2 f+°° 



fc2 

Voir), (56) 
' ~ (D-l p'' 



$(r) 



dpJoip) e 



(57) 



F(r) 



and ((57)) can be integrated numerically. In this latter, 
Jo{p) = sinc(p) = sm{p)/p is the Bessel function. For 
small values of the radial coordinate "r" (large values of 
"p" ) we get 



$ w -2GnM (const.) r^''+^, 



(58) 



where const, sa 3 x 10^ tt for A = 1 and Gn = 1- The 
potential ([58)) is regular for r — > and a plot of the exact 
potential (|57)) for 7 = 3 and M = 10 is given in FigHJ 





FIG. 1: Plot of the gravitational potential in D = 4 for 7 = 3 
and M = 10 (A = Gjv = 1). 



Making use of the gravitational potential (|57)) . for the 
sake of simplicity we can study the "Newtonian cosmol- 
ogy" in Z? = 4. To derive the Friedmann equation in 
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Newtonian cosmology, we need the kinetic and potential 
energy of a t est p article and we must implement energy 
conservation [5lL Is^. We now consider an observer in 
a uniform expanding medium of mass density p. Be- 
cause the Universe is homogeneous and isotropic, we can 
assume any point to be its center. We then identify a 
particle of mass m at a radial distance r. Due to New- 
ton's theorem, the particle only feels a force from the 
material at smaller values of r. The material has mass 
M = Airpr^ /3 and the constant total energy of the test 
particle is E = T+U, where T = mP /2, U = m<I>(r) and 



$(r) is given in ([57)1 . Because the Universe is homoge- 
neous, we apply this argument to any couple of particles, 
which allows us to introduce comoving coordinates de- 
fined by r = a{t)x. For the same reason, the real distance 
r is related to the comoving distance x by a(t), which is 
a function of time alone. When dividing E = T + U hy 
a{t)'^x'^, we get the modified Friedmann equation j53l[54 



H2 



— ^pF(a)--, (59) 



where K = -2E/mx^ and F(a) is defined in dST]). To 
maintain homogeneity, the quantity E must depend on 
the comoving coordinates according to E oc x'^. For the 
same reason, we have rescaled (x 1/x^ obtaining an 
equation independent from x so that homogeneity is re- 
spected. Hereon we assume K = 0. As wc know the scal- 
ing of the gravitational potential for small values of "r" , 
we then can get the Friedmann equation near a{t) w 0. 
In this limit we find 



H 



2 SttGn 2(const.) 27+2 

3 TT 



(60) 



where the constant is defined right after ((58)) . It is clear 
that this Universe is singularity-free, since H — > when 
a{t) — 0. Furthermore, if we add the cosmological con- 
stant Acc to ([59]) . then the equation (|60)) reads 



(61) 



^2 _ SttGn 2 (const.) ^27+2 _j_ 



This shows that the cosmological constant dominates the 
Universe at high energy whatever kind of matter we in- 
troduce, either dust or radiation. Consequently, the Uni- 
verse follows a natural de Sitter evolution at the Planck 
scale. We can numerically integrate equation ([59)) as 
FigUl shows for both dust matter and radiation. 

In a ZJ-dimensional spacetime, the Newtonian cosmol- 
ogy can be theorized in a way similar to the 4-dimensional 
case. The Friedmann equation (for K = 0) together with 
the fluid one, reads 





FIG. 2: Plot of H^(a) in D = 4 for 7 = 3 and A = Gat = 1 
in a universe dominated by dust matter on the left and radi- 
ation on the right. The dashed lines account for the classical 
trajectories. 



dE + P dV = T dS [P is the pressure, T the temperature 
and S the entropy) and assuming reversible expansion 
{dS ~ 0). In D dimensions, the state equation for dust 
matter is P = and for radiation is p ^ {D — 1)P, so 
that the fluid equation implies 



Pdust 



Fdust 



Prad. 



Prad. 



(63) 



When a ^ 0, the function F£,(a) « {aKf^+'^^+'^ and 
the Friedmann equation simplifies to 



H2 



N 



{D-2){D-l) 



p{aA) 



2y+2N+2 



(64) 



In (|M)) goes to zero for both radiation and dust mat- 
ter, which means that the D-dimensional Newtonian cos- 
mos are singularity-free. 

Consistently with the singularity-free cosmology we 
have illustrated so far, the black hole solutions turn out 
to be regular, as we are going to demonstrate in Z) = 4. 

Following [55l - [6^ . the equations of motion for the 
above theory (up to square curvature terms) are 



G 



0{V^Rpa) = 8nGNV{z)Tf,,, (65) 



where the argument z ~ —Da as defined throughout in 
the paper^. 

Since we are going to solve the Einstein equations ne- 
glecting curvature square terms, then we have to impose 
the conservation {V {z)Tfj_^) = in order for the theory 
to be compatible with the Bianchi identities. Conversely, 



167rG 



N 



{D-2){D-l 



■pFc(a), 



p+{D-l)^-j{p + P)^Q (fluid cq.), (62) 

with ¥o{a) as defined in ([55)1 . The fiuid equation in (|62l) 
has been obtained using the first law of thermodynamics 



In general, a differential equation with an infinity number of 
derivative has not a well-defined initial value problem and it 
needs an infinite number of initial conditions. It is shown in 
|66|468|| that in a general framework each pole of the propagator 
contributes two initial data to the final solution. This is pre- 
cisely our case because the only pole in the bare propagator is 
the massless graviton and the theory has a well defined Cauchy 
problem. 
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the exact equations of motion satisfy the Bianchi identi- 
ties because the theory presents general covariance. The 
condition V^{V{z)T^ij) = compensates the truncation 
in the modified Einstein equations (|65p. 

Our main purpose is to solve the field equations by 
assuming a static source, which means that the four- 
velocity field has only a non- vanishing time-like com- 
ponent u^" = (u",0), u° = [79l - [98t . We consider 
the component o of the energy-momentum tensor for 
a static source of mass M ([53|) . In polar coordinates, 
TOq = p = M(5(r)/47rr^ The metric of our space- 
time is assumed to be given by the usual static, spheri- 
cally symmetric Schwarzschild form 



ds^ = F{r)dt^ - 



dr^ 



F(r) = 1 - 



Fir) 
2GNm{r) 



(66) 



The effective energy density and pressures are defined by 



V{z)Tt' 



Diag(p°,-P;,-Pl,-Pl). (67) 



SttGtv 

We temporarily adopt free-falling Cartesian-like coordi- 
nates [6J, |65| to calculate the effective energy density, 
assuming pj+N+i{z) = in (|38l). 



p%x) 



v{-a^)T% 

d^k 



MVi-DA)d{x) 



(68) 



= M 



(27r)3 
2M '■+°° 



(2^) 



-H( 



p sin(p) dp, 



where r — \x\ is the radial coordinate. Here we intro- 
duced the Fourier-transform for the Dirac delta function 
and we also introduced a new dimcnsionless variable in 
the momentum space, p = \k\r, where k is the physical 
momentum. The energy density distribution defined in 
(|68|) respects spherical symmetry. We evaluated numeri- 
cally the integral in and the resulting energy density 
is plotted in FiglS] In the low energy limit we can ex- 
pand H{z) for z = — D/A^ <C 1 and we can integrate 
analytically (|68| 



p%r) 



2M 



(27r)2r3 



+00 



p sin(p) dp. (69) 



The result is extremely complex, and its plot is given in 
FiglSj however, the Taylor expansion near r « gener- 
ates a constant leading order p'^(r) oc MA^. 

The covariant conservation and the additional condi- 
tion, = —grr fully specify the form of V{z)T^i, and 
the Einstein's equations reads 



dm[r) 

dr 
1 dF 
F~dr ' 
dPl _ 
dr 



2Gn (m(r) 



AnP^r^) 



r(r 
1 dF 
'2F~dF 



-2GNm{r)) ' 

(p^ + P,^) + -(PI-P,°). (70) 





FIG. 3: Plot of the energy density for m = 10 in Planck units 
assuming A — mp. The solid line is a plot of (|68|) without any 
approximation. The dashed line refers to the energy density 
profile ([SSI) in the limit -D/A^ < 1. 



Because of the complicated energy density profile, this is 
how the first Einstein equation would fit in (|70|) 



m(r) = 47r 



p^(r'). 



(71) 



However, the energy density goes to zero at infinity, re- 
producing the asymptotic Schwarzschild spacetime with 
m{r) sa M (constant). On the other hand, it is easy to 
calculate the energy density profile close to r w since 
H{z) -> log 2^ for z +0O (or r in ([HH])). In this 
regime, m{r) oc M A^r^ and, for a more general mono- 



mial p^-|_i(z) = z 



7 + 1 



i{r) oc M(Ar)2')'+2. The function 



F{r) in (|66p near to r w is approximated by 

P(r) « 1 - c Gat M A^^+^ r^^+\ (72) 

where c is a dimcnsionless constant. 

We show now that the metric has at least two horizons: 
an event horizon and a Cauchy horizon. The metric in- 
terpolates two asymptotic flat regions, one at infinity and 
the other in r = 0, so that we can write the g~j^ = —F 
component in the following way 



Fir) 



1 



2MG(r) 



(73) 



Here G(r) Gn for r -> oo, G(r) cx Gn r'^''+'^ for r -> 
and G(r) does not depend on the mass M. The function 
F{r) goes to "1" in both limits (for r +oo and r 
0). Since M is a multiplicative constant, we can always 
vary it for a fixed value of the radial coordinate r, such 
that F{r) becomes negative. From this it follows that 
the function F{r) must change sign at least twice. The 
second equation in ([70)1 is solved by = —p^ and the 
third one defines the transversal pressure once the energy 
density p'^ is known. For the lapse function F{r) in ([7^. 
we can calculate the Ricci scalar and the Kretschmann 
invariant 



R = cGn M A2^+2 (27 + 2) (27 + 3) r^''-!, (74) 
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By evaluating the above curvature tensors at the origin 
one finds that they are finite for 7 > 1/2 and, in par- 
ticular, for the minimal supcr-rcnormalizable theory in 
D = 4 with 7^3. 




FIG. 4: The first plot shows the function F{r) for the energy 
profile (|69|) and H{z) defined in H39|) with the parameter 7 — 
3. The ADM mass values are M = 1 and M = 10 (in Planck 
units) for the first and the second plot respectively. 



The form factor V{z) is able to tame the curvature 
singularity of the Schwarzschild solution at least for the 
truncation of the theory here analyzed. However, we be- 
lieve that the higher order corrections to the Einstein's 
equations will not substantially change the fundamental 
remarkable feature of the solutions found in this section 

Besides, we can exactly (but only numerically) inte- 
grate the modified Einstein equations (p5)) for the en- 
ergy density defined in (|68p . Using the integral form of 
the mass function (j7ip , we achieve the metric component 
F{r) defined in ((66|) . The numerical results are plotted in 
FigHlfor different values of the ADM mass M. The met- 
ric function F{r) can intersect no times, twice or more 
than twice the horizontal axis according to the value of 
the ADM mass M. This may enable "multi- horizon black 
holes" as an exact solution of the modified equations of 
motion ([55)1 . 

We expect the same features to be maintained in any 
dimension D > A. 



We can now address a more general class of theories 
following the Efimov's book on non local interactions [l^] • 
Let us consider the propagator in the following general 
form 



zA2 



(75) 



(the notation is rather compatible with the graviton 
propagator ([M)) and z — Da)- 

As was shown by Efimov , the nonlocal field theory 
is "unitary" and "microcausal" provided that the follow- 
ing properties are satisfied by V{z), 



I. V{z) is an entire analytic function in the complex 
z-plane and it has a finite order of growth 1/2 ^ 
p < +00 i.e. 36 > 0, c > so that 



\V{z)\ ^ ce''!"!' 



(76) 



IL When Re(z) — >• -|-oo (fc^ — > -l-oo), V{z) decreases 
quite rapidly. For example, we can consider the 
following cases. 

a. V{z) = O ( |Jpr) , a > For a = ^ the 

theory is not super-renormalizable, but may 
still be renormalizable. 



\z\^\V{z)\ = 0, ViV > 0. 



b. limRe(z)->+o 

m. [v{z)]* = v{z*). 

IV. V{0) = 1. 

V. The function V{z) can be non-negative on the real 
axis, i.e. V{x) > 0. 

Here are some examples of possible functions: 

A. Va{z) = e^^" for ?i e N+, 
the weight is p = ?i < +00, 



B. Vb(.) = (^)' 



C. Vc{z) = rr{l + s)^ {s>0). 

When V{z) = Vb(z) or Vc{z), the functions hi{z) in the 
action are not entire functions, so they do no longer meet 
our minimal requirement. 

A more refined growth measure is obtained by defining 
the order jo(0i, ^2) for V{z) in the angle 9i ^ argz < ^2- 
It is a remarkable property of entire functions that, for 
appropriate V{z), p{0i,02) may range from zero to ar- 
bitrarily large values as 9i,92 change. The function 
V{z)^^ = expiJ(z) that we have introduced and exten- 
sively studied in the first part of the paper exhibits at 
most polynomial behavior along the real axis and it is of 
infinite order p ~ +00 in the full complex plane. 

To expand on the point II. a, we calculate the propa- 
gator in the coordinate space for a general form factor 
V{z). The Fourier transform of ([75]) reads 



G(x) 



(2^) 



D 



£ee1/A, (77) 



where we neglected any tensorial structure and we as- 
sumed Euclidean signature. Changing the existing coor- 
dinates into D-dimensional spherical ones and integrat- 
ing ((77|) in the angular variables, we get 



G{x) 



D-3 
TT 2 



(2^) 



D-l 



du 



u—viue) 



D - 1 



(78) 
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where we have introduced the variable u = k'^ . From II. a, 
V{ul'^) = 0(l/u") for u — >■ +00 and since o-^i ~ const, 
for 0, the propagator in the coincidence limit is 

finite only for certain values of a, 



G(0) oc 



+00 



g-4 

duu 2 



< 00 



a > 



D-2 



(79) 



For _D = 4 the two-point function in the coordinate space 
is 



G(x) 



1 



i2ny 



-\-oo 



duV{ul')^l^^, (80) 
2v ux'^ 



where " Ji" is the Bessel function of the first kind "J„(z)" . 
Using II. a {u — >■ +00) and/or the short distances limit 
— >■ 0, the propagator ([50]) reads 



O 



Gix) 



0{\n{x^)) for a 
0(1) for a > 1 . 



for < a < 1 , 
1, 



(81) 



Only for a > 1, G(0) < +00 in the coincidence limit. This 
is further evidence that super-renormalizability requires 
a > {D — 2)/2. Later on we will show that we may still 
have renormalizability for a = (D — 2)/2 in the case study 
D = 4. 

We now move to investigate the general theory (|15p 
with form factor Va{z) — exp(— z"), which satisfies the 
property II. b., in the entire functions hi{z) defined in 
(Plj) . The high energy propagator reads 



(82) 



■ Dr^h - 



The m-graviton interaction has the same scaling, since it 
can be written in the following schematic way 

h . . . , (83) 

where = rj'^^d^dij. The notation ". . ." indicates other 
sub-leading interaction terms coming from the covariant 
D'Alembertian operator. Placing an upper bound to the 
amplitude with L-loops, we find 




I-V 



(84) 



In the last step we used again the topological identity 
/ = V + L — 1. The L-loops amplitude is UV finite 



for 1/ > 1 and it diverges as "fc^" for L = 1. Only 1- 
loop divergences survive in this theory. Therefore, the 
theory is super-renormalizable and unitary, as well as 
microcausal as pointed out in (27l - l32j . 

To calculate the gravitational potential for n = 1, 
it suffices to replace exp H {p^ / A^) exp(p^/r^A^) 
within the integral (IMl) . The result is: 



V(0 = TT 



T 



kM 
2 47r^ 
D-3 



1 



E„ 



-T 



D-2, r"^ 



(85) 



To prove the regularity of the graviton solution, we ex- 
pand ([55]) near r = 0, so that we get the following finite 
leading term 



V(o) 



{D-?,)iT— 
For £> = 4, (155|) simplifies to 
GnM 



$(r) 



Er 



r A 



(86) 



(87) 



The gravitational potential is regular in r = and its 
value is $(0) = -GnMK/ ^Jtt. For rt > 1, the potential 
is still regular in r = and it takes the value $(0) oc 
—GjyMA with a slightly different coefficient. 

In the case n = 1, wc can always solve the equations of 
motion (j65p for a spherically symmetric _D-dimensional 
spacetime with metric 



dsl) 



FDir)dt' 



dr^ 



Fair) 



r dVlD-2: 



where d^D-2 is described in terms of Z? — 2 angles. The 
form factor V{z) — cxp — z gives a smearing of the source 
and the energy density reads 



p-{r)^V{z)T\^M 



A2 
47r 



(89) 



Integrating the "00" component of the modified Einstein 
equations (|55|). we get the function Fair), 



Foir) 
'D 

7 ■ 



2MGd 



7 



1 r2A2 



2 ' 4 
T{{D-l)/2) = 



^^^^4 dt o . 





D - 1 



t 



-t~ 



- 1 



! for D odd , 



{D-'i)\\ 

2iD-2)/2 



for D even. (90) 



where [Gd] = AP^^ ■ The other components in ((65)) are 
solved by p'^ = —P°, while the covariant conservation of 
the effective energy tensor V{z)T^i, determines 



{D - 2)-i r drp" , i^l,...,D-2. (91) 
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The metric has a "de Sitter core" near the origin r 
where 



AMGd a 



D-l 



{D - l)2^-27r(-D-3)/2 



(92) 



from which descends a singularity-free spacetime. All 
the other properties of the metric have been extensively 
studied in 

The repercussions of this study affect several fields as 
emerges from pr evious inv estigations in LHC black hole 
phenomenology [lOOl . llOlj . gauge _gravity duality jl02l | 
and early universe cosmology |62l . |97| . Specifically it has 
been shown that the resulting blac k hole tend to emit 
softer particles on the brane |lOl| . a fact which is in 
marked contrast with previous results based on classical 
metrics. 



Another special theory we wish to explore is defined 
by the following form factor, 



V{z)-. 
H{z) 



1 



(93) 



[7i. + r(0,z2^+2)]+log[z^+^ 



Re(z2^^+2) >0. 

This form factor has been achieved from ()38p by choosing 
7 = 0. The theory satisfies all the properties I - V of the 
second class of theories here examined. In particular, the 
behavior of the entire functions hi{z) for \z\ — > +oo is. 



lim |/t,(z)| ^ for z in: (94) 

|z|— >+oo 

C = {z I - 9 < argz < +9 , tt - 9 < argz < tt + 9}, 
for 9 = 7r/(4A^ + 4). 

Since in even dimensions N — (-Devon ~ 4)/2, the entire 
functions hi{z) in _D = 4 approache a constant for |z| — > 
+00. This theory embodies the quadratic Stelle action 
in the ultraviolet limit but without any ghost pole in the 
propagator. The form factor cross-connects the quadratic 
action in the infrared with an equivalent theory in the 
ultraviolet. 

The theory in question meets the property II. a for the 
critical value a = (D - 2)/2 (a = 1 in D = 4) The 
amplitudes are divergent at each order in the loop ex- 
pansion and the maximal superficial degree of divergence 
from ([7]) or © is 5 = D as it occurs in the local theory. 
Therefore, the theory ceases to be super-renormalizable, 
but it preserves renormalizability and unitarity as it can 
be inferred from (j46p and (j47p with the entire function 
H{z) defined in 



The gravitational potential can be obtained integrating 
([57)) with the form factor (^5)) . The potential is regular 
everywhere and 4>(r) — 14 GatMA^ r near the point r = 
0. Because the metric scaling ^(7-) « 1 — (const.) MA^ r 
in r = 0, black hole solutions are not singula rity-free, as 
proved by the diverging curvature invariants [I05l4l07j . 

Let us assume that the coupling constants in D = A 
satisfy the following relation (see ([50]) *) 



bo{Zbg - 1) = 3ao(l - Zao) + 1 



(95) 



in the ultraviolet regime. Then the Lagra ngian turns out 
to be conformal invariant at high energy [104| . 



3 



(96) 



However, the same result can not be achieved for the lo- 
cal Stelle's theory, because the relationship between the 
coupling constants, for which we get a conformal invari- 
ant action, is the same one by which the theory loses its 
renormalizability. 




FIG. 5: Plot of 1^2(2)^ with the form factor defined in (|93} 
for D = A and then A*" = 0. To draw this plot we have taken 
= 2, A = 1 and feo = 1 in (ED. 



A fundamental quantity that explains "the spacetime 
structure" is the spectral dimension (hereafter dg). This 
is not only a tool to compare different approaches to 
quantum gravity, but it is actually a device to extract 
information about the physics within the spacetime. It 
is equivalent knowing either the spectral dimension, the 
propagator or the gravitational potential. Schematically, 



Gravitational Potential. 



(97) 



See page 147 and pages 246-252 of the Efimov's study [27 



More explicitly, if we know the spectral dimension, then 
we also know the heat kernel (see below), from which we 
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can derive the propagator and the ensuing gravitational 
potential. Clearly, the reverse relationship is true as well. 

Here we calculate the spacetime spectral dimension 
flow from short to long distances for the three different 
cases already discussed in the paper, 



Form Factor 1. Vi{z) 



Form Factor 2. ¥2(2) = e^""" , n G N+ , 

Form Factor 3. V,iz) = eM7.+r(o,.--)]+iog(."^^)^ 

As we arc going to show, renormalizability, along with 
unitarity, implies a spectral dimension dg < 1 for the 
the form factor 1, ds = for the form factor 2 and 
ds = 2 for the form factor 3. Let us recall the definition 
of spectral dimension in quantum gravity. Such defini- 
tion is borrowed from the theory of diffusion processes 
on fractals and adapted to the quantum gravity con- 
text. In the Brownian motion of a test particle mov- 
ing on a I?-dimensional Riemannian manifold A4 with a 
fixed smooth metric g'^^^ix), the probability density for 
the particle to diffuse from x' to x during the fictitious 
time T is the heat-kernel K{x,x';T). This satisfies the 
heat equation 



dTK{x,x';T) 



N$K{x,x'-T), 



(98) 



where A®? denotes the usual covariant Laplacian at low 
energy, which may undergo substantial modifications in 
the ultra-violet regime. In particular, we are interested 
in the effective Laplacian at high energy on the flat back- 
ground ((7°y = 77^1/) where the graviton propagates. The 
heat- kernel is a matrix element of the operator exp(T Ag) 
acting on the real Hilbert space of position eigenstates 



K{x,x'-T) = {x'\cMT^f)\x). 



(99) 



Its trace per volume unit. 



PIT) 



_ Jd"xy/^K{x,x;T) 



V 



Tr exp(TA5?) 



V 



(100) 



can be interpreted as an average return probability. Here, 
V = J d^x denotes the total volume. It is acknowl- 
edged that P{T) possesses an asymptotic expansion for 
OoftheformP(T) = (47r T)-^/^ A„ T". The 

coefficients An have a geometric meaning, i.e. Aq is the 
volume of the manifold. Knowing P{T), one can recover 
the dimensionality of the manifold M. as the limit for 
large T of 



dlnT • 



(101) 



This formula defines the fractal dimension we are going 
to use. 

Omitting the tensorial structure in p4p . which does 
not affect the spectral dimension, we can easily obtain 



the heat-kernel. We know that the propagator (in the 
coordinate space) and the heat-kernel are related by [toI - 
M 



r+00 

Gix,x') = / dTK{x,x';T) 
Jo 



(102) 



\2t^) Jo 



where 



G{x,x') 



d^k 



is the Fourier transform of ([M)) . By inverting (|102p with 
respect to the heat-kernel in the momentum space, we 
get 



i^(fc;T)cxe-^'='^('='/^')" 



(103) 



The necessary trace to calculate the average return prob- 
ability is obtained from the Fourier transform of (jl03l) , 

Kix,x';T)o: [ d^fc e'^'^' ^('^'Z^')" e''^^^"^'). (104) 



Now we are ready to calculate the average return proba- 
bility defined in (|100p as follows 



P{T) cx 



(105) 



We then proceed to calculate explicitly the spectral di- 
mension for the three different form factors listed above. 

* Form factor 1. At high energy, V{k^)-^ ~ ^2^+2^+2. 
Therefore, we can calculate the integral (jlOSp and then 
the spectral dimension defined in (|10ip for small T is 



P{T) cx T 



^D/{2-i+2N+i) 



D 



7 + iV + 2 



(106) 



Since the parameter 7 > Devon/ 2 or 7 > (fodd — l)/2, 
the spectral dimension is < 1 V Z? and 



lim dg 

D-^+oo 



1 



(107) 



applies, which is a "universal" property of this class of 
theories. Using the explicit form of the entire function 
H{k^ /A^) given in (|40)) . we calculate the spectral dimen- 
sion at all energy scales as the fictitious time T varies. 
Integrating numerically (|105p . we can plot directly the 
spectral dimension achieving the graphical result in Figl6] 
for D = 4, 6, 8, 10 and 7 = 3,4, 5, 6, respectively. 

* Form Factor 2. For simplicity we consider the case 
n = 1, even though for n > 1 the result is qualitatively 
the same. Given the form factor in the momentum space 
V^(fc2/A2) = exp(— /s^/A^), the propagator scaling reads 



o{kr' 



(X 



(108) 
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FIG. 6: Plot of the spectral dimension as a function of the 
fictitious time T for D = 4, 6, 8, 10 and the minimal values 
7 = 3, 4, 5, 6 in (|40p . The graph on the right shows that the 
spectral dimension approximates ds = 1 when increasing D. 





FIG. 8: Plot of the spectral dimension as a function of the 
fictitious time T for D = 4, 6, 8, 10 and form factor Vz{z). 



spectral dimension in the ultraviolet regime. In other 
words, once perturbative renormalizability and unitarity 
have set the form factors, the spectral dimension in the 
ultraviolet regime flows to the same "critical point" or 
"accumulation point" independently from the topological 
dimension D. From this evidence, we can infer that any 
consistent theory of quantum gravity must satisfy the 
following fractal property. 



ds ^ 2 VZ) in the ultraviolet regime. 



(Ill) 



FIG. 7: Plot of the spectral dimension as a function of the 
fictitious time T for D = 4,6,8, 10 and form factor V{z) = 
exp{—z). 



and the heat-kernel can be calculated analytically. 



K{x, x'; T) = 



g 4(T+1/A-!) 

[4^ (r + A-2)]^ 



(109) 



as verifiable by (|108|) . Applying (|101|) . we find that the 
spectral dimension is 



ds = 



T 



D. 



(110) 



which clearly goes to zero for T — > and approaches 
da = D for T — > +CX). A plot of the spectral dimension 
flow is given in Figl7]for D = 4, 6, 8, 10. 

* Form factor 3. In this case the spectral dimension in 
the ultraviolet regime is = 2, VZ?. A plot of the spec- 
tral dimension flow is given in Fig |8] for Z) = 4, 6, 8, 10. 

It is remarkable to note that for all the three classes 
of theories parametrized by (7, n, D) and studied in this 
section, we always flnd an accumulation point for the 



CONCLUDING REMARKS 

This study is a synthesis of concepts coming from non- 
local quantum field theory [l^], p article ph ysics, general 
relativity and string field theory |l23l4l28| . In this ar- 
ticle we suggested ways to quantize gravity, relying on 
the perturbative approach that has been so successful 
for the other fundamental forces. We introduced a non- 
local extension of the higher-derivative gravity, which is 
perturbatively renormalizable and unitary in any dimen- 
sion D. The four-dimensional theory is easily obtained 
from the Stelle theory [13] by introducing in the action 
two entire functions, a.k.a. "form factors" , between the 
Ricci tensor square and the Ricci scalar square. 



Rho{-D/A^)R, 
R^, /i2(-n/A2)i?'^^ 



(112) 



In the multidimensional spacetime we preserved the 
two "delocalization-operators" as in (|112p and we im- 
plemented a finite number of local operators required 
(and/or generated) by the quantum consistency of the 
theory. These local operators 02n{dg) contain 2n- 
derivatives of the metric tensor up to the mass dimension 
[02n{dg)] ^ M^. The action may also present other ir- 
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relevant operators, whose couplings constants have neg- 
ative mass dimension. The full action reads 



(113) 



N 



-J2['^nR (-Qa)" R + bn i?^. (-Da)" R^' 

n=0 

-Rhni-DA) R + i?^.. /i2(-nA) i?^" 



0{R^ 



R 



N+2 



Finite number of terms 



0{R''+-') + 0{R{-UKy'+' R) 
^ . ' 

Non rcnormalizcd operators 



The main reason for introducing the entire functions 
h2{z) and hQ{z) is to avoid ghosts (or rather the 
poltergeists: states of negative norm) and any other new 
pole in the graviton propagator. The unitarity require- 
ment implies the following entire functions. 



h2{z) 



ho{z) = - 



2 ^ 

-o(i?(-nA)^+ii?)), 



o(i?(-nA)^+ii?)). 



(114) 



The first set of operators in the last line of (|113p is sub- 
ject to renormalization at quantum level, whereas the 
second set remains classical as it can be proved by power- 
counting arguments. Clearly, the non-renormalizcd op- 
erators 0{R (— □a)^"''^ i?)) in (|114p can be eliminated in 
both the action and the entire functions. 

The form factors V{z)~^ studied in this paper can es- 
sentially show two possible high energy behaviors, either 
polynomial or exponential. In the first case, the operators 
0{R^^^) may affect the renormalizability of the theory, 
therefore the polynomial asymptotic degree of V{z)~^ 
has to be increased. In the second case, the same lo- 
cal operators do not thwart the renormalizability of the 
theory at all. 

Let us gather here all the quantities to be measured to 
define the theory, 

(a„,6„), OsCjis^iV, (115) 
0(c„ i?") , 3 ^ n ^ + 2 (local relevant operators). 
Viz) , non-local form factor , 

0{dn i?") , ri iV + 3 (local irrelevant operators). 

As pointed out throughout the paper, the irrelevant op- 
erators can always be introduced in any physical theory 
as long as they do not invalidate its unitarity and renor- 
malizability. The most important step we need to take 
is to assess whether the physical measurable quantities 
are affected or not by such irrelevant operators in (|115|) . 
If we can assume that no physical quantity is susceptible 
to such operators, we can then empirically infer that the 
coupling constants equal "zero" . 



The question lingers whether the form factor V{z) is 
measurable or not. In principle we can treat Viz) as 
one of the form factors for the scattering of the nucleus 
by electrons. In a gravitational theory, such measure 
represents the graviton scattering amplitude, as well as 
the modifications to the gravitational potential or the 
light-bending. The four-gravitons amplitude will have 
the general structure 



Ag. = -44g. (s, t, u; F(s, i, u); £1,2,3,4) 



(116) 



where €1,2,3,4 are the four gravitons polarizations and 
s,i, u the Mandelstam variables. Since V^z) has to be 
an entire function, we can falsify the theory by compar- 
ing the experimental four-gravitons amplitude with the 
theoretical prediction (|116p . A final remark about the 
tree-level unitarity of the theory has to be put forward at 
this point. Although we consider this issue still open 0- 
0,11^, we must also acknowledge that, at high energy, the 
total cross section of a nonlocal interacted theory must 
not exceed that of a local one. Intuitively, the reason is 
that "nonlocal particles" must manifest the transparency 
property at high energy because of their non-zero size. 
For example, if the amplitude in the momentum space 
grows exponentially 



A 



(117) 



the total cross-section will satisfy the following upper 
bound fH, 



a'tot(s) sS const, s'' ^ log^(s) , 



(118) 



which grows logarithmically in the ultraviolet regime if 
p=l. 

In this paper, we studied the following three classes of 
form factors V{z), 

Form Factor 1. Vi{z) = e^^^^^^ , 

Form Factor 2. ^2(2:) = e"^" , n e N+ , 

Form Factor 3. V3(z) = e~^"'=o(^) , 

where 

H,{z) = \ [7s + r(0,z2'^+2A.+2)] +log(,7+A^+l). (119) 

We systematically showed the power-counting renormal- 
izability and the tree-level unitarity. The theories de- 
fined by the form factors Vi{z) and V2{z) result to be 
renormalizable at one loop and finite from two loops up- 
ward. More precisely, the theories turn out to be super- 
renormalizable because the covariant counter-terms have 
less derivatives then the classical action and the coeffi- 
cients of the terms with more derivatives do not need 
any kind of infinity renormalization as synthesized in the 
first part of this section. However, we argu e that a su- 
persymmetric extension of the theory [122i | can make it 
finite at one loop as well. For the third choice V^{z), 
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the theory is merely renormalizable, and no other pole 
beyond the graviton one appears in the propagator. 

We solved the linearized equations of motion and we 
proved that the gravitation potential is regular in r = 
for all the choices of form factors compatible with renor- 
malizability and unitarity. We also included Black hole 
spherical symmetric solutions omitting higher curvature 
corrections to the equation of motions. For two out of 
three form factors {Vi{z) and V2{z)) the solutions are 
regular and the classical singularity is replaced by a "de 
Sitter-like core" in r = 0. For the third choice V3(z), 
black holes are still singular, although the divergence is 
attenuated. 

For Vi{z), we proved that the "Newtonian cosmology" 
is singularity-free in any dimension D and the Universe 
spontaneously follows a dc Sitter evolution at the "Planck 
scale" for any matter content (either dust or radiation), 
since the cosmological constant dominates the effective 
energy tensor at high energy. In a £)-dimcnsional space- 
time the modified Friedmann equation (for K = and 
Acc = 0) reads 



N 



{D-2){D-1) 



(120) 



where Foia) w {aAf'>+^^+'^ for a w 0. 

Finally, we have provided an extensive analysis of the 
spectral dimension for any D and for the three classes of 
theories. In the ultraviolet regime, the spectral dimen- 
sion takes on different values for the three cases: 



Vi{z) 



ds < 1 



4 = 0, 



where H^{z) is defined in (|119|) . Once the class of theories 
compatible with rcnormalizability and unitarity is de- 
fined, the spectral dimension has the same short-distance 
"critical value" or "accumulation point" for any value of 
the topological dimension D. This is a "universal" prop- 
erty of the theories here studied. 

We would like to conclude this section by identify- 
ing some similarities between the second class of super- 
renormalizable theories and "string field t heory". U sing 
the results found at the e nd of the Eighties jl23l - ll28j and 



several more recent ideas [491, Il29l | , the string field theory 
has the following schematic structure for the spacetime 
bosonic and fcrmionic fields, 



S 



d^x ( i(/),A'y(n)<^, - iHjk4>;Mk ] , (121) 



where 



, ln(3 



%/3 , 



Kij (□) « □ for open as well as close bosonic strings, and 
a' is the inverse mass square in string theory. By a field 
redefinition |l29{ . the action (|12ip simplifies to 



(122) 



We can immediately observe that the kinetic term in 
(|122p has the same scaling of the linearized theory stud- 
ied in this paper for the exponential form factor V2{z) = 
exp(— z) (n = 1). If we expand (|113p in powers of the 
graviton field neglecting the exponential factor in the in- 
teraction, the thrcc-graviton vertex is quite similar to the 
one in (|122|) . However, the general covariancc in (|113p 
implies the same leading scaling in the kinetic term as 
well as in the interaction vertexes and we are unable to 
get a finite theory at any order in the loop expansion. As 
already pointed out in this section, one possible loophole 
to this puzzle could be a supersymmctric extension of the 
action in ((T5|) . 

About the finiteness of string theory, we are likely to 
endorse the following ideas. Due to the presence of the 
exponential factor, the effective string theory in (|122p 
manifests an asymmetry between the kinetic and the in- 
teraction terms. Contrary to our covariant action (jll3p . 
such asymmetrical state implies that the string theory 
does not manifest any divergence. The well-known "soft- 
ness" of the high energy tree-level amplitudes also de- 
scends from the same asymmetry. 

However, the comparison here proposed can only be 
qualitative and partial because, unlike the effective string 
field theory, ours is a general covariant theory. Indeed, 
general coordinate invariance in string theory can only 
be achieved through cancellations amo ng co ntributions 
from infinitely many interactions terms |126| . However, 
we do not exclude that a supersymmetric extension of our 
theory (jllSp can be framed within "M-theory" as one of 
its possible vacuums. 



APPENDIX: 3D HIGHER-DERIVATIVE 
QUANTUM GRAVITY 

In this section, as a particular toy model, we consider a 
nonlocal gener aliza tion of the SZ^ higher derivative grav- 
ity studied in |109| . The nonlocal action is 



S 



g\ [aR + R a{n)R + i?^,/3(n)ii'^n , (123) 



where the two "form factors" «(□) and /?(□) are "entire 
functions" of the covariant D'Alembertian operator. We 
introduce the following definitions, 

a{U) := ao + h^i-Ui,) , /?(□) := /3o + h^i-U!,) , (124) 

where Da := D/A^ and ft-Oi ^2 arc entire functions. The 
two form factors have dimensions: [Q;(n)] = [/3(n)] = 
LP'. The Lagrangian, complete with the gauge fixing and 
ghost terms, reads 



£ = £<, + £ 



GF 



C. 



GH, 



(125) 
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where Cg is expressed by (|123|) . and the graviton fluctu- 
ation h^^ is defined by 



(126) 



Imposing the BRST invariance on the full Lagrangian 
(|125p , we ean get the gauge- fixing and ghost terms of the 
action. The BRST transformation for the fields in ()125p 
appears 



5bc^ = 



(127) 



where c^tC^, are the anti-commuting ghosts fields, Bp 
is the auxiliary field, 5\ is an anti-commuting con- 
stant parameter and uj{U\) is a weight function of the 
d'Alembertian operator. The dimensions of the fields 
are: [h] = [c] = [c] = L'^'-i/^, [B] = £"3/2, 

[6\] = L^^y . The BRST transformation for the graviton 
field in (|126p can be extracted from 

= SKr'dpC'' + r'dpc'' - -g^-dpcP - dp~g>^- 

= 6\V^^''cP, (128) 

which implies Ssh'"' ^ kSWi^" cP [io|. The gauge fix- 
ing and ghost Lagrangian can both be expressed as a 
BRST variation of the following functional 



£gf + Cgh = iSs [cp.{dyhP- 
= -BpLo{U)d^ht"' 



aBP/2)]- (129) 



-B^oj{U)B^. 



To obtain the graviton propagator, we first eliminate the 
auxiliary field i?^ to get the following gauge fixing La- 
grangian, 



and then we assemble the quadratic part of (|125p , namely 

(131) 



+ncj(n) p(iV« + P*"'''((8a(n) + a /?(□)) - a) 

_^2p(o^-')((8 «(□) -f 3 /?(□))□ - (T + w(n)/a) + \/2 X 



In (|13ip we have introduced the 31) projectors. Using 
the orthogonality and the completeness property of the 
projectors, we find the graviton propagator 



D{k) 



(27r)3 fc2 



p(2) 



+ 



(8Q(fc2) 3/3(fc2))fc2 



3(1) 



P(0.-) V2 



In the harmonic gauge d^h^^ = (or a = 0), the propa- 
gator considerably simplifies to 



D{k) 



1 



(27r)3fc2 



p(2) 

Ml) 



p(0,s) 



(132) 



where the following notation has been introduced, 

h2{z) -.^ p{z)zA^ ~ a , 

ho{z) (8a(z) + 3/3{z))zA^ + a, 

z -Da. (133) 

As in the D-dimensional case ([HT]) . we choose the entire 
functions h2{z), ho{z) compatibly with the properties (i)- 
(iii), 



h2{z) 



a{Viz)-^-l)+PoA^ 
A^z 



M»)^ ^'"^'""-''-^°'°"" . (134) 

8A^z 

Assuming the theory to be renormalized at a particular 
scale /io; we identify 



ao = ao(Mo) , /3o = /3o(mo)- 



(135) 



In this case h2 — ho ~ V{z) ^ and the propagator sim- 
plifies to 



D{k) 



V 



(27r)3 a fc2 



[p(2) _ p(0,s)j 



(136) 



The pole structure of the propagator is the same one as 
that in the local theory, because V{z) is an entire function 
with no zeros in the complex plane. In Z? = 3 there are 
no local degrees of freedom, and therefore the amplitude 
(Uni) is identical to zero. 

What we have presented in this section is a toy model 
of modified nonlocal gravity that we think it might be 
interesting to expand on in connection with other three- 
dimensional theories studied in recent years |llCll - ll21 1 . 
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